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The phase diagram of QCD with two and four flavors: results with HYP 
fermions * 

A. Hasenfratz a and F. Knechtli a t 

a Department of Physics, University of Colorado, Campus Box 390, Boulder, CO 80309, USA 

We study the finite temperature phase transition of four and two flavor staggered fermions with hypercubic 
smeared link actions. These preliminary studies suggest that the improved flavor symmetry of the fermionic 
action can have a significant effect on the finite temperature phase diagram. 



1. INTRODUCTION 

In this talk we present preliminary results of 
the two and four flavor QCD finite temperature 
phase diagram on Nt — 4 temporal lattices us- 
ing a hypercubic blocked (HYP)Q,j^] staggered 
fcrmion action. Because of the large lattice spac- 
ing these simulations are more relevant to illus- 
trate the effect of the HYP action than to extract 
physical values. Our goal here is to show how the 
improved flavor symmetry of the HYP fermions 
influences the phase diagram. 

In any fermionic action the thin link gauge con- 
nections can be replaced by some sort of smeared 
links. Smearing partially removes short range 
vacuum fluctuations (dislocations) improving the 
topological properties of the gauge field configu- 
ration and improving the chiral properties of the 
fermionic action. At the same time smearing can 
distort the short range properties of the config- 
uration introducing new type of lattice artifacts. 
An optimal smearing is a transformation where 
short range distortions are minimal while most 
of the dislocations are removed. The hypercubic 
blocking was designed to be very compact, and 
its parameters are optimized non-perturbatively 
to remove most of the vacuum fluctuations. In- 
terestingly the non-perturbatively optimized pa- 
rameters remove the perturbative tree level flavor 
symmetry violating terms of the action. 
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The HYP blocking can be combined with any 
fermionic action. Here we consider staggered 
fermions. In [[l] we showed in the quenched 
approximation that flavor symmetry with HYP 
staggered fermions is improved by about an or- 
der of magnitude relative to thin link staggered 
fermions. The improved flavor symmetry is espe- 
cially important for physical quantities that are 
sensitive to chiral symmetry. The finite temper- 
ature phase transition restores chiral symmetry 
at vanishing quark mass. Theoretical predictions 
!| for this transition are based on the full chi- 
ral and flavor symmetry of the fermionic action. 
Lattice simulations studying the nature of the fi- 
nite temperature phase transition at small quark 
masses could therefore be very sensitive to chiral 
symmetry violations of the lattice action. Since 
the finite temperature phase transition occurs at 
large lattice spacing, a w 0.3 — 0.1 5fm on typical 
temporal lattices Nt = 4 — 8, the use of chirally 
improved actions could be essential there. 

Before presenting our numerical results we 
mention a puzzle about the flavor dependence 
of the finite temperature phase transition. Phc- 
nomenological instanton model calculations pre- 
dict a strong flavor dependence for the chiral tran- 
sition temperatures [B. For Nf = 2 flavors the 
chiral transition is expected to be second order 
at T c » 150MeV with crossover at finite quark 
mass values. For Nf — 3 flavors the transition 
is predicted to be first order, T c m lOOMeV in 
the chiral limit while for Nf — 4 flavors the chiral 
transition, if it exists at all, is at a very small tern- 
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perature. At finite quark mass both the Nf = 3 
and 4 flavor cases show first order phase transi- 
tion that persists up to some critical quark mass 
value and the phase transition temperature de- 
pends strongly on the quark mass. Lattice sim- 
ulations observe only a weak flavor dependence, 
especially for Nf = 4 flavors. Instanton models 
rely on plausible but unproven assumptions yet 
they predict the zero temperature chiral behav- 
ior of QCD very successfully It is puzzling why 
their prediction is different from lattice calcula- 
tions in the finite temperature case. 

2. N F = 4 FLAVOR SIMULATIONS 

We performed simulations on Nt = 4 tempo- 
ral lattices with N$ = 8, 10 and 16 spatial size. 
The details of the simulations can be found in 
!| . To see the effect of smearing first we compare 
thin link, one level of APE and HYP smeared ac- 
tions at approximately matched quark masses. In 
quenched simulations we found that the three ac- 
tions are matched at quark masses amwYP = 0.1, 
oroAPEi = 0.08, and amthin = 0.06. Dynamical 
simulations show that once the quark masses are 
matched, the lattice spacing at fixed gauge cou- 
pling P remains approximately matched also. As 
it can be seen in figure [l] the strong first order 
transition observed with thin link action weakens 
and disappears as smearing is introduced at con- 
stant physical quark mass. One might wonder if 
the HYP action might show a first order tran- 
sition at some smaller gauge coupling. We can- 
not exclude this possibility without actual sim- 
ulations, but even if there is a first order phase 
transition, it occurs at a much smaller temper- 
ature than the phase transition of the thin link 
action suggesting that it is not a physically rele- 
vant transition. 

We did not find first order signals even at 
smaller quark masses close to the physical light 
quark mass indicating that our simulations were 
all above the end point of the first order transition 
line. To identify the crossover it is traditional to 
look at the susceptibility of the quark condensate. 
We found only very broad peaks in the suscepti- 
bility. However the peak of the susceptibility does 
not necessarily identify the crossover transition. 
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Figure 1. The chiral condensate < ipip > for 
the thin link, APEl and HYP actions at approx- 
imately matched quark masses (Nf = 4). 

It is more sensitive to nearby critical points like 
the end point of the first order line. If the tran- 
sition line has a strong quark mass-temperature 
dependence, this could be very different from the 
actual crossover region. Because of that concern 
we decided to use the chiral condensate at fixed 
gauge coupling to identify the crossover temper- 
ature. At small quark masses we expect the con- 
densate to depend linearly on the quark mass, 

< iptb > £ 

= h C 

m q m q 

where £ and c are constants independent of the 
quark mass. We distinguish three different cases. 
If all quark masses are in the chirally symmet- 
ric phase, X = 0. If all quark masses are in the 
chirally broken phase, £ > is the chiral conden- 
sate. If at fixed gauge coupling the smaller quark 
mass points are in the chirally symmetric while 
the larger quark mass points are in the chirally 
broken phase, the formers predict £ = while the 
latter ones predict a non-zero condensate at the 
transition point. £ however is not this value but 
its linear extrapolation to zero quark mass and 
could take any value, even negative ones. Figure 
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Figure 2. < ?/)^ > /am q as the function of \/am q 
at three different gauge coupling values (iV/ = 4) . 
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Figure 3. The predicted crossover line for the 
Nj = 4 flavor system on Nt — 4 temporal lat- 
tices. 

H illustrates such an extrapolation at j3 = 4.8, 4.9 
and 5.0. The (3 — 5.0 points are all in the chi- 
rally symmetric phase, the other two are mixed. 
From the plot we can read off the transition mass 
at fixed gauge coupling. The negative slope at 
larger quark masses indicates that the condensate 
at the crossover is small. The predicted crossover 
region is plotted in figure |. At = 5.0 the 
crossover occurs above am q = 0.1 corresponding 
to a temperature around 220MeV. At j3 = 4.95, 
arrit = 0.04 the crossover temperature is approx- 
imately T — 170MeV. Whether this temperature 
changes significantly as we move to smaller quark 
masses has to be investigated in the future. In 
any case our data indicate that the first order end 
point occurs at a fairly small quark mass value 
(and possibly small temperatures). 
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Figure 4. < rptp > /am q as the function of the 
gauge coupling for Nf = 2 flavors at two different 
quark mass values. 

3. Nf — 2 FLAVOR SIMULATIONS 

We have preliminary results of the two flavor 
system on 8 3 x 4 lattices at several quark masses. 
In the simulations we approximate the square 
root of the fermionic determinant with a finite 
polynomial. Some of the details of the algorithm 
are discussed in || . The phase diagram does not 
show significant deviation from the thin link ac- 
tion results. Figure H shows < xptp > /am q as the 
function of the gauge coupling for am q = 0.01 
and am q = 0.04. The figure suggests a phase 
transition around (3 = 5.2. Further data at quark 
mass am q = 0.02, 0.03 and 0.06 support this con- 
clusion. The phase transition for the two flavor 
system on Nt = 4 is a crossover, it occurs around 
(3 = 5.2 and depends only weekly on the quark 
mass. The susceptibility of the chiral condensate 
predicts the same value. 
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